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Л.В. ЧЕРКЕСОВА 


АНАЛИЗ РЕШЕНИЙ НЕЛИНЕЙНОГО УРАВНЕНИЯ МАТЬЕ-ХИЛЛА 
С УЧЁТОМ СОБСТВЕННОГО ТОКА РЕЗОНАНСНОГО КОНТУРА 


Дан анализ решений нелинейного дифференциального уравнения Матье-Хилла, описывающего нелинейный 
резонатор или нелинейно-параметрическую зонную (или пазонную) систему (НПС), работающую в высших 
зонах неустойчивости колебаний на высших гармониках. Получены устойчивое и неустойчивое решения это- 
го уравнения с учётом влияния собственного тока резонансного контура НПС на модуляцию его индуктивно- 
сти. 

Ключевые слова. резонансный контур, нелинейный резонатор, нелинейно-параметрическая зонная систе- 
ма, амплитуда колебаний. 


Введение. Нелинейный параметрический эффект широко используется в электронике для создания 
малошумящих усилителей и высокоэффективных генераторов, элементов памяти, логических элемен- 
тов — регистров, счетчиков, дешифраторов, сумматоров и функциональных преобразователей и др. 
Применение электронной аппаратуры, имеющей в своем составе нелинейные резонаторы, чрезвы- 
чайно широко. Нелинейные резонаторы представляют собой колебательные параметрические зонные 
(или пазонные) системы, работающие в первой или высших зонах неустойчивости колебаний, на 
высших гармониках. 

Резонансные схемы, использующие колебательные контуры, представляют собой нелиней- 
но-параметрические зонные (пазонные) системы (НПС). Они широко используются в качестве вы- 
сококачественных усилителей, генераторов, преобразователей, а в цифровой технике -— тригге- 
ров, логических элементов и элементов памяти. Их применение позволяет не только улучшить 
технические характеристики электронной аппаратуры, но и повышает их надежность в условиях 
неблагоприятного воздействия окружающей среды. Достоинства таких систем обусловливают их 
использование в измерительной технике, в качестве датчиков НЧ-диапазона как малошумящих функ- 
циональных устройств. 

Постановка задачи. Изучение свойств и закономерностей явлений в НПС, работающих в высших 
зонах неустойчивости, весьма актуально и возможно при разработке точных методик и широкого 
класса экспериментальных методов. Подобные исследования лежат в плоскости наноэлектроники при 
изучении нелинейных параметрических резонансных явлений в наноразмерных структурах. 

Методы и результаты исследования. Нелинейно-параметрические взаимодействия распро- 
странены в природе очень широко и присущи многим нелинейным средам. Колебательные про- 
цессы в НПС удобно моделировать с помощью нелинейных электрических цепей, в которых легко 
получить, описать и исследовать все явления и процессы, связанные с нелинейно- 
параметрическими взаимодействиями. Нелинейная пазонная система индуктивного типа в общем 
случае описывается нелинейным дифференциальным уравнением [1]: 








А? | . 
оба, да.) х ид =0, й 
а, (г) 4 Ц Ко) 
ба и К ЕН Ш а 1 Е 
где 9(,(г))= В) + О) и (КП) КО) ы ес соответственно 
коэффициент затухания и квадрат частоты собственных колебаний резонансного 


контура. 
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Для решения уравнения (1) предположим, что влияние собственного тока резонансного 
контура параметрического преобразователя на рассматриваемую систему отсутствует. Тогда оно 
преобразуется к уравнению: 








4? Аи) т ь 
=. +20(0) +9(х К =0. 
Используя подстановку 1(г) = Ух = се ‚ его можно преобразовать к виду: 
йо 
РОЖИ =0, те РЗ -6°-° Ю 


При исследовании НПС индуктивного типа кривую намагничивания сердечников катушек 
индуктивности часто считают безгистеризисной и аппроксимируют гиперболическим синусом [1, 
2]. Тогда ток в цепи накачки в случае, если параметрический преобразователь подключен к гене- 


ратору с напряжением (0 (г) =И„, соз@{, определяется следующим образом: 
| [971 фо 
1, (=—х5й ЕЙ 

И’ 25” 


В работах [1, 2] было показано, что изменение индуктивности одного сердечника резо- 
нансного контура будет иметь вид: 








у 
= Е. 
сИ(Ь $ ( @г)) 
10 | ВОт 
где [о =——— - индуктивность резонансного контура при КЁ=0; Б= — приведен- 
[7] 25’@& 


ная амплитуда тока накачки; —И/- число витков обмотки накачки; $, [ — площадь 
поперечного сечения и длина средней линии сердечника; о — круговая частота на- 
пряжения накачки. 

На практике большой интерес представляет НПС индуктивного типа, в которой затухание 
очень мало, т.е. сопротивление активных потерь значительно меньше реактивных сопротивлений: 
(А<<м, К<< №). 

Предположим, что К = 0 (потери в системе о В этом случае: 


Г/Р). [60а _ ВЮ 





(= е Е? (в) = = 05 хей(Ьзш а), 3 

ПО 10’ оС т. } 8) 
1 

где 5 ме: представляет собой квадрат собственной частоты резонансного контура 
0 


НПС при отсутствии модуляции индуктивности. 
Функция Е? (1) имеет физический смысл квадрата круговой частоты собственных коле- 


баний резонансного контура. В дальнейшем для решения уравнения (2) с учётом (3) используем 
следующее разложение: 


сй(Ьзш 8) = [3 (Ь) + у (-1)* 1. (6)хсоз(2060), 


= 
где ББ) - модифицированные функции Бесселя чётных порядков. 
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В случае К = 0 уравнение (2) имеет вид: 
“и › и 
= + (1(6)+2» (1 х 1, (Б)соз(2КеЮ)И = 0. 
К= 
Сделав замену переменной о{=т, получим уравнение Хилла в виде: 


ти + [*) х | (6) + У (-1/* 15, (6) комкор =0. (4) 
@ 





: 
Чт = 


Для решения уравнения (4) можно применить метод Уиттекера, согласно которому 
[е.0) 
необходимо, чтобы ряд 15% (5) <® сходился [3]. 
к=0 
Докажем сходимость этого ряда. Известно, что 


(2) + у (—-1)* Г., (2) = с08(2), 


где [/(2) - модифицированные функции Бесселя порядка п [4]. 
Прибавив в последней формуле к обеим частям [(2), получим: 


2». (—-1)* 1. (2) = с08(2) + (2). 


Заменив в последнем выражении 2 на 1х, а также используя известные формулы 


1, @х) =РТ, (х) и соз(&х) = сй(х), окончательно получаем: 


2 г: сй(х)+ 1 (х) 
2169 Е ЗЕ . 


[.0) 
Таким образом, ряд УР 2к(Ь) <со сходится. 
К=0 
Применим метод Уиттекера для решения уравнения (4). В случае малости коэффици- 
ентов при сомножителе со$(2Кт) можно получить достаточно хорошее приближение, исполь- 
зуя три центральных строки и три центральных столбца определителя Хилла бесконечного 
порядка [3]. 
И ‘0 
Для уравнения Хилла Е (т)хИ =0 чаще всего рассматриваются два случая, в 


одном из которых функция Е” (г) имеет вид: 


Е*(т)=0, +20 соз(2т) +20, соз(4т)+..., (5) 
[е.0) 
где коэффициенты 0, - заданные постоянные, иряд »`@„ абсолютно сходится. 
п=0 


Положив 0..=6„ , ищем решение уравнения Хилла в виде: 


[.9) 
т 2тт 
а лаем. 
й=—с 
где и — характеристический показатель. 
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Для определения амплитуды а, имеем систему уравнений: 


[©.6) 
(и+2т)?а, + Убалт=0 (и=....-Ъ 0,1... (6) 
т=—ёо 
Для определения ги составляется определитель Хилла (1/4). Известно [5], что корни опре- 
делителя Хилла являются корнями уравнения: 


р 2 
и? [| - А) ЕЙ -- | 


(7) 


поэтому необходимо найти (0), а затем из (7) найти д. 


Для количественного анализа решения уравнения Хилла ограничимся только коэффи- 
циентами ду, 01, 02, так как при большом количестве коэффициентов 0„ решение становится 
неоправданно громоздким при незначительном возрастании точности. Отбрасывание остатка 
ряда (5) при вычислениях производится исходя из условия 0 >> 09т+1 (Т21). 

Рассмотрим случаи, когда приведённая амплитуда тока накачки принимает следующие 
значения: В =1; 1,6; 2,2; при о. = - 

о 2 

В случае В =1, о =2, находим следующие значения: (0) =1,002; и = -0,563/. Подставив 
найденное значение д в (6) и положив 0.; = 0;, получим: 

—4х(п- 0,282}? ха, +6, х (а +а:)+6, ха, =0. 

Приняв значения л равными 0, 1, —1, получим однородную систему из трёх уравнений с пятью 
неизвестными. Положим а› = а.2= 0. Тогда получим однородную систему трех уравнений с тремя неиз- 
вестными, в которой свободной переменной будет а. Решив эту систему, получим: 

а1=- 0,005 а, а =- 0,019 ао. 

Таким образом, величину амплитуды колебаний определить нельзя. Это подтверждает тот 
факт, что теория линейных дифференциальных уравнений с периодическими коэффициентами не 
может дать ответа на вопросы о величине стационарной амплитуды, процессе ее установления и 
т.д. [1, 2, 6]. 

Положим аи = —1. Тогда: а.,= 0,005, а, = 0,019. Подставив значения ах, ат, а, и взяв мнимую 
часть, находим искомую функцию Ц (т): 

((г)= зщ (0,563: )+ 0,019 х зи (1,4375) 0,005 х зй1(2,563т). 


Ток (т) равен: 


п-1 


сИ(Бзшт) 


( 
С учетом (8) окончательно получаем выражения для тока (т) во всех случаях. 
При Б = 1и о = 20 получим: 


ЦТ) = О(т)х (8) 


Кг) = —@.269 $11(0,563т) + 0,161х $11(1,437т) — 
0 
— 0,143 х 511 (2,5631) — 0,005 х $11(3,437т) + 0,003 х 11(4,563т)). 
При Б = 16ио = 26%: 
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Кг) = — 0,78 511(0,67т) + 0,52 х $11(1,337) — 
0 
— 0,39 х51(2,67т)-0,05х 511(3,33т) + 0,02 х з11(4,67т)). 
При В = 2,2 и о = 260: 


Кг) = — (2,97 х и1(0,87т) + 1,8 1х (1137) —$1(2,87т) — 


0 
— 0,39 хз1п(3,13т) + 0х з11(4,87т) + 0,03 жзш(513т)). 


На рис.1 приведены в нормированных единицах графики функций 1=Кт) для значений 
Ь=1, 16и 2,2. Из графиков видно, что в областях устойчивого решения уравнения Хилла 
установившиеся параметрические колебания существовать не могут. В этих областях на- 
блюдаются случайные колебания с постоянно изменяющейся амплитудой и фазой, т.е., как 
говорят, в этих областях нелинейный резонатор НПС индуктивного типа «шумит». 





Рис.1. Кривые, описывающие механизм усиления тока в резонансном контуре при увеличении 
приведенной амплитуды тока накачки Ь. Обозначения на рис.: Б=1 - сплошная линия; 
Ь=1,6 - пунктирная линия; Ь=2,2 — штрихпунктирная линия 


В работах [1, 2] был рассмотрен случай питания нелинейно-параметрического преобразовате- 
ля индуктивного типа от генератора напряжения. Аналогично производится расчет и в случае пита- 
ния от генератора тока. Рассмотрим случай, когда генератор имеет синусоидальную форму тока: 


йе (е) = 1, $1(00). (9) 


165 


Физико-математические науки 








Динамическая индуктивность контура описывается выражением [1, 2, 7]: 
257? 1 1 


‚ ШН)=—х ===. (10) 
ВУИ” +а7Р а пм: 
[9 


Подставив выражение (9) в формулу (10) и сделав преобразования, получим выражение 
для динамической индуктивности: 





(= 











2Г 
Ки= 212 - | (11) 
ИТ 
у ибо 
а 
2 
5 
где 7 — ов! ` 
212 
Обозначим 5 ы =2т2. Применив соотношение 2 $1 * (69) =1-с05$(2627) и сделав 
Га 


несложные преобразования, получим: 


—> 


7. 
ЕЙ 12 
| _ т? с03(24) (12) 
1+7? 
20 
1+ т? 
В случае А=0 (для системы без потерь) для квадрата круговой частоты собственных колеба- 


ний резонансного контура К * (+) выше была получена формула (3). Тогда с учётом выражения 
(12): 


Г = 


где 9 = 








1 2 
Е? (р) = 5 +т' }- т = с0$(201). 
т 


Отсюда, в случае питания параметрического преобразователя от генератора тока, урав- 
нение (4) преобразуется к виду: 


2 2 2 

ее 5% 1+ т? 1 - Е 

ат 21 о 1+т 

Воспользуемся известной формулой биноминального ряда: 

п(п-1) и п(п-1(п-2) о 
! 3! 

Применив эту формулу и подставив значение для переменной х 











(1+х)" =1+их+ 








_ т? с0$(25) 
1+4 
получим следующее уравнение для определения напряжения (/: 
2 
РИ [а р 1 т? с0о$(2т) 1[т? с0з(25) 
+—| — | №14 | 1-—. ...|ХО =0, 1 
ат? :“ 2 шт 8 1417 м 
В уравнении (13) ряд, представленный в квадратных скобках, сходится для любых значе- 
й т? соз(2т) 
ний 77 и т так как |-————^5—^ <1. 
1+т 
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@) х 
Рассмотрим решение уравнения (13) для значений: —*=0,9; [, =103.10°А (при 
@) 


т=2,444). После подстановки исходных данных получим: 
40 


2 
т 





+ (1,0204 —0,458 со5(2т) —0,046 соз(4т)-...)хИ = 0. 


Ограничимся для анализа следующим уравнением, в котором для Е? (1) учтены только два 


слагаемых ряда, так как большее количество слагаемых приводит к громоздким вычислениям при 
незначительном увеличении точности: 
а’ 
= (1,0204 - 0,458 с05$(2т))хИ =0. (14) 
ат 
Найдем неустойчивое решение уравнения (14), соответствующее первой области не- 
устойчивости, используя метод Уиттекера [3], согласно которому общее решение уравнения 
Матье 
2 
Чу 
2 
т 





+ (а-24с05(2т))ху=0 


может быть записано в виде: 
у=Ае"" хФ(т,0)+ Ае ""хФ(т,-ов)), 
где А, › — произвольные постоянные; д - характеристический показатель; о- параметр, прини- 
И 
мающий для неустойчивых решений значения между0и ^ 5 | 
Функция Ф(т, о) имеет для первой области неустойчивости следующий вид: 


Ф(т, о) = зш(т-о)+ 53 зщ(3т 0) + С3 соз(3т -©) +... 


3 : 
где 5, -- 59+ 49° 008(27)+.., а. 49 (20) +... : 


Для определения с имеем уравнение: 


Я ое аи Ба +1-а=0. (15) 
4 64 64 1536 8 
Для уравнения (14) 4=0,229. Из уравнения (15) находим, что о =-0,872. Тогда для д, $3, 
Сз получим следующие значения: /=0,113; $ =-0,029; С; =-0,002. Ограничимся точностью вы- 
числений, равной Зх 10° (считаем С; =0), тогда: 
Ф(т, с) = зш(т + 0,872) 0,029 зв (Зт + 0,872). 
Отсюда находим функцию у(?: 


(г) = Але [ит + 0,872) — 0,029 зи(3т + 0,872) |+ 


0,113 [< . и 
+ де 13° [те - 0,872) — 0,029 зи(3т — 0,872) | 
Считая, что у(0) = 0, из (16) получаем А, = 45. 
А 45 ы - 
Постоянные — ‚—= имеют физический смысл амплитуды начальных колебаний в ре- 
о ^0 


зонансном контуре НПС индуктивного типа. Для реальных систем, работающих в первой об- 
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е А А -. ь т 
ласти неустойчивости, можно принять = Тогда, учитывая, что АЕ. 
Ро т) 


получаем: 





2 
2 
Кг) = а Пе ЕЕ. те т). 
21 1+т 


Взяв два слагаемых ряда (при 17 = 2,444), окончательно получим для тока Кт): 


Кт) = - х ео [оС + 0,872) —0,029 (3х + 0,872)|+ 
и (17) 
че 13" [а(т — 0,872) — 0,029 5(3т — 0,872); (10,428 с0$(2т))х10-°. 

Полученное решение (17) не учитывает влияние собственного тока резонансного контура 
НПС на модуляцию индуктивности [(6) этого же контура. Это не позволяет математически точно 
описать процессы установления колебаний в резонансном контуре системы после достижения ими 
некоторой величины, определяемой, в основном, нелинейностью параметрической системы. 

Для устранения этой проблемы воспользуемся методом последовательных приближений, опи- 
санным в [8, 9], и сучетом влияния собственного тока резонансного контура параметрической систе- 
мы на модуляцию его индуктивности, получим выражение для суммарной индуктивности: 

Цт,Кт)) = у х 258 
[, +КоР? О В 





Величина 27 с учетом первой гармоники (как имеющей наибольшую амплитуду по срав- 
нению с ВЫСШИМИ гармониками) будет иметь вид: 
2 
й АА (е"" +е #7) 
ре 20 


2_ И? 
Ра? 
В последних полученных выражениях ток Их) может как суммироваться, так и вычитаться 
с током /[». Это обусловлено схемными особенностями параметрических преобразователей индук- 
тивного типа, у которых резонансная обмотка на одном сердечнике включена согласно с обмот- 
кой накачки, а на другом встречно [1, 2, 10, 11]. Учитывая это, суммарная индуктивность /Г(т, (т)) 
и величина 2? будут равны: 


1 





2т (18) 


И? 1 1 
+ 
В | н+коРь +Ро? ры -кеРи2 + Ра? 
ры Е 
Ра?" р 


Це, е))= 








(19) 





2т ср? ит 


Коэффициент 4 в (18) определяется в начале каждого периода из выражения 


Я = 1+ [(0)Р ‚ где т(О) - величина т, вычисленная в точке, где /Г„ = 0. Такое вычисление 


коэффициента @ вполне оправдано, так как этот коэффициент в течение полупериода изменения 
индуктивности изменяется медленно. 
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Вычисление (т) по выражению (17) с учётом влияния собственного тока резонансного 
контура системы производится до момента перехода решения уравнения (13) из области неустой- 
чивости в область устойчивости. Дальнейшее вычисление (т) производится по устойчивому ре- 
шению уравнения (13). 

На рис.2 и рис.3 приведены расчетная и экспериментальная кривые нарастания и уста- 
новления колебаний тока в резонансном контуре параметрического преобразователя индуктив- 
ного типа, имеющего следующие конструктивные характеристики: 30 витков в обмотках накач- 
ки и резонансных обмотках; сердечники К7 х 4х 2 из феррита 1500 НМ и ёмкость конденсатора 
резонансного контура С = 0,2 мкФ. 


250 


-250 





Рис. 2. Расчетная кривая нарастания тока в резонансном контуре 





Рис. 3. Экспериментальная кривая нарастания тока в резонансном контуре 


Установившееся значение тока резонансного контура для данного параметрического пре- 
образователя рассчитывалось по выражению, являющемуся устойчивым решением уравнения 
(13) с учетом (19): 


Кт) = ты (5а() — 0,026 зи(3т))х (1- 0,428 с05(2т)), (20) 
2 
где 1 - амплитудное значение тока собственных колебаний контура в момент перехода решения 
из области неустойчивости в область устойчивости. 
Аналогично рассмотренному выводу выражения (17) при питании параметрического пре- 
образователя от генератора тока для случая его питания от генератора напряжения такое реше- 
ние имеет вид: 


Кг) = Е фо [1(е + 0,898) — 0,014 зп(3т + 0,898) |+ 
и (21) 
+е°'55* х[зш(е — 0,898) —0,014 (35 — 0,898) (1,26-0,28 соз(2т))х10-3. 
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Выводы. Таким образом, при нахождении только устойчивого решения уравнения Матье-Хилла 
методом Уиттекера амплитуду колебаний в резонансном контуре НПС индуктивного типа опреде- 
лить не представляется возможным. Для её определения необходимо найти и учесть неустойчи- 
вое решение уравнения. Сопоставляя выражения (17) и (19) для тока при питании параметриче- 
ского преобразователя от генератора тока и генератора напряжения, следует отметить, что ток 
в резонансном контуре нелинейно-параметрической зонной системы нарастает быстрее в случае 
её питания от генератора тока. Выражения (20) и (21) определяют решение нелинейного диффе- 
ренциального уравнения Матье-Хилла с учётом установившегося значения собственного тока ре- 
зонансного контура для изучаемого параметрического преобразователя. 

Результаты работы представляют значительный теоретический и практический интерес 
для проектирования радиоэлектронных систем. 
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